AEDC-TR-75-2 


^tB  0  ^ 

r7  f| 


=^®333333I> 


REAL  GAS  SCALE  EFFECTS  ON  HYPERSONIC  LAMINAR 
BOUNDARY-LAYER  PARAMETERS  INCLUDING 
EFFECTS  OF  ENTROPY-LAYER  SWALLOWING 


VON  KARMAN  GAS  DYNAMICS  FACILITY 
ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


December  1975 


Final  Report  for  Period  March  21  —  October  20,  1974 


pT'rerty  of  U.  S.  A:‘r  Foro@ 

AElJ-'O 

F4-0600-75-C-0001 


Approved  for  public  release;  distribution  unlimited. 


Prepared  for 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION  (JSC) 
JOHNSON  SPACE  CENTER 
HOUSTON,  TEXAS  77058 


NOTICES 


When  U.  S.  Government  drawings  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement 
operation,  the  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,'  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  pr  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense 
Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be 
considered  in  any  sense  as  an  endorsement  of  the  product  by  the  United 
States  Air  Force  or  tlie  Government. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 


APPROVAL  STATEMENT 


This  technical  report  has  been  reviewed  and  is  approved  for  publication 


FOR  THE  COMMANDER 


JIMMY  W.  MULLINS 
Lt  Colonel,  USAF 

Chief  Air  Force  Test  Director,  VKF 
Directorate  of  Test 


FRANK  J.  PASSARELLO 
Colonel,  USAF 
Director  of  Test 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

AEDC--TR-75-2 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  ("and  SubddeJ  5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

REAL  GAS  SCALE  EFFECTS  ON  HYPERSONIC  Final  Report,  Mar  21  ~ 

LAMINAR  BOUNDARY-LAYER  PARAMETERS  INCLUD-  Oct  20,  1974 
ING  EFFECTS  OF  ENTROPY-LAYER  SWALLOWING 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORfs; 


8.  CONTRACT  OR  GRANT  NUMBERfs; 


J.  C.  Adams,  Jr.,  W.  R.  Martlndale, 

A.  W.  Mayne,  Jr.,  and  E.  O.  Marchand, 
ARO,  Inc. 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

.  V  AREA  a  WORK  UNIT  NUMBERS 

Arnold  Engineering  Development  Center  (XO) 

Arnold  Air  Force  Station,  TN  37389  Program  Element  921E 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  12.  REPORT  DATE 

National  Aeronautics  and  Space  December  1975 

Administration  (NASA/JSC) , 

Houston,  TX  77058 


14.  MONITORING  AGENCY  NAME  a  ADDRESSfi/ d/Herenl  Irani  Conirolffnfi  Otiice)  15.  SECURITY  CLASS,  (ol  this  report) 


13.  number  of  PAGES 

50 


UNCLASSIFIED 


1S«,  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

N/A 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  In  Block  20,  if  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


Available  in  DDC. 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 


space  shuttles 
laminar  boundary  layer 
nonequilibrium  flow 
equilibrium  flow 


boundary-layer  transition 
nose  bluntness  effects 
aerodynamic  heating 
hypersonic  flow 


20.  abstract  ^ConUnue  on  reverse  side  if  necessary  and  Identify  by  block  number) 

Inviscid  and  viscous  (laminar  boundary  layer)  flow-field  calcu¬ 
lations  under  perfect  gas  hypersonic  wind  tunnel  and  equilibrium 
real  gas  flight  conditions  are  presented  for  the  windward  center- 
line  of  the  Rockwell  International  139  Space  Shuttle  Orbiter  at  30 
deg  angle  of  attack.  Correlation  parameters  for  laminar  boundary- 
layer  edge  quantities  and  surface  heat  transfer  are  developed  which 
properly  account  for  entropy-layer-swallowing  effects  under  both 


FORM 
1  JAN  73 


EDITION  OF  1  NOV  65  IS  OBSOLETE 


UHCLASSIFIED 


UNCLASSIFIED 


19,  Continued 

hypersonic  reentry  flight 
viscous  hypersonic  flow  fields 


20,  Continued 

real  and  perfect  gas  conditions.  A  cursory  examination  of 
chemical  nonequilibrium  effects  on  the  inviscid  flow  field 
based  on  oblique  shock  relaxation  is  presented.  Some  im¬ 
plications  of  the  proposed  correlation  parameters  on 
boundary-layer  transition  are  discussed. 


AFS  c 

Arnold  AFS  Tenn 


UNCLASSIFIED 


AEDC-TR-75-2 


PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at 
the  request  of  the  National  Aeronautics  and  Space  Administration  (NASA/ 
JSC),  for  the  Rockwell  International  Space  Division.  The  results  pre¬ 
sented  were  obtained  by  ARC,  Inc.  (a  subsidiary  of  Sverdrup  &  Parcel 
and  Associates,  Inc. )  contract  operator  of  AEDC,  AFSC,  Arnold  Air 
Force  Station,  Tennessee,  under  ARO  Project  Numbers  VC523  and 
V43Y-09A.  The  authors  of  this  report  were  J.  C.  Adams,  Jr.  ,  W.  R. 
Martindale,  A.  W.  Mayne,  Jr.  ,  and  E.  O.  Marchand,  ARO,  Inc. 

The  analysis  was  completed  on  October  20,  1974,  and  the  manuscript 
(ARO  Control  No.  ARO-VKF-TR-74- 104)  was  submitted  for  publication 
on  October  25,  1974. 

Acknowledgement  and  appreciation  are  extended  to  D.  E.  Schramm 
of  Rockwell  International  Space  Division  for  providing  the  flight  tra¬ 
jectory  and  wall  temperature  information  used  in  this  study. 


1 


AEDC-TR-75-2 


CONTENTS 

Page 

1.0  INTRODUCTION .  7 

2.0  FLOW -FIELD  MODEL .  8 

3.0  FLOW-FIELD  CALCULATION  METHODS 

3.  1  Boundary  Layer  and  Inviscid  Flow .  10 

3.  2  Chemical  Nonequilibrium  Streamtube  .......  14 

4.0  FREE-STREAM  CONDITIONS 

4.1  Wind  Tunnel .  17 

4.2  Flight .  17 

5.0  FLOW-FIELD  MODEL  VERIFICATION  UNDER 

HYPERSONIC  WIND  TUNNEL  CONDITIONS  .......  20 

6.0  CORRELATION  OF  BOUNDARY-LAYER  PARAMETERS 

INCLUDING  EFFECTS  OF  ENTROPY-LAYER 

SWALLOWING .  27 

7.0  ESTIMATE  OF  CHEMICAL  NONEQUILIBRIUM 

EFFECTS  ON  OUTER  EDGE  PARAMETERS .  36 

8.0  IMPLICATIONS  RELATIVE  TO  BOUNDARY- 

LAYER  TRANSITION . 42 

9.0  CONCLUDING  SUMMARY .  43 

REFERENCES . 44 

ILLUSTRATIONS 

Figure 

1.  Rockwell  International  139  Space  Shuttle  Orbiter 

Configuration . 8 

2.  Comparison  Between  139  Orbiter  Windward  Centerline 

Surface  Profile  at  30  deg  Angle  of  Attack  and  a  31 -deg 
Asymptotic  Half-Angle  Hyperboloid . 9 

3.  Hyperboloid  Geometry,  Showing  Inviscid  Streamline 

Swallowing . 12 

4.  Chemical  Nonequilibrium  Oblique  Shock  Relaxation  ...  15 

5.  Rockwell  International  SSV  Trajectory  No.  14040  ....  19 

6.  Comparison  Between  Hyperboloid  Surface  Pressure 
and  Shock  Shape  Calculations  and  AEDC/VKF  Tunnel 

B  f^ata  ...................0.0.  21 


3 


AEDC-TR-75-2 


Figure  Page 


7.  Comparison  Between  Hyperboloid  Boundary -Layer 
Edge  Parameter  Calculations  and  AEDC/VKF 

Tunnel  B  Data  ......  .o  ...  . .  22 

8.  Comparison  Between  Hyperboloid  Boundary-Layer 

Profile  Calculations  and  AEDC/VKF  Tunnel  B  Data  ...  23 


9.  Comparison  Between  Hyperboloid  Surface  Heat- 
Transfer  Calculations  and  AEDC/VKF  Tunnel  B 
13ata  ........  .......  .........  26 

10.  Correspondence  of  Physical  Surface  Distance,  s,  for 

a  Given  Value  of  the  Similarity  Parameter,  S .  28 

11.  Logarithmic  Plot  of  the  Laminar  Boundary-Layer 

Outer  Edge  Mach  Number  Correlation  Parameter  ....  30 

12.  Logarithmic  Plot  of  the  Laminar  Boundary -Layer 
Outer  Edge  Unit  Reynolds  Number  Correlation 

!Paraxn.eter  ...........  ..........  31 

13.  Logarithmic  Plot  of  the  Laminar  Boundary -Layer 

Thickness  Correlation  Parameter  . .  33 

14.  Logarithmic  Plot  of  the  Laminar  Boundary-Layer 

Surface  Heat-Transfer  Correlation  Parameter .  36 

15.  Chemical  Nonequilibrium  Oblique  Shock  Relaxation 

Effects  on  the  Local  Mach  Number .  38 

16.  Chemical  Nonequilibrium  Oblique  Shock  Relaxation 

Effects  on  the  Local  Unit  Reynolds  Number .  40 

TABLES 


1.  COORDINATE  CORRESPONDENCE .  10 

2.  CHEMICAL  REACTIONS . 16 

3.  REACTION  RATES . 16 

4.  THIRD-BODY  EFFICIENCIES  . .  16 

5.  TUNNEL  B  CONDITIONS  . .  17 


4 


AEDC-TR-75-2 


Page 

6.  MODEL  WALL  TEMPERATURE  RATIO  DISTRIBUTIONS 

UNDER  TUNNEL  B  CONDITIONS .  18 

7.  ROCKWELL  INTERNATIONAL  SSV  TRAJECTORY 

NOv  14040  .  18 

8.  CENTERLINE  RADIATIVE  EQUILIBRIUM  WALL 

TEMPERATURE  DISTRIBUTIONS  FOR  ROCKWELL 
INTERNATIONAL  SSV  TRAJECTORY  NO.  14040  ...  .  20 

9.  CONICAL  MACH  NUMBER  BASED  ON  INVISCID 
SOLUTION  FOR  A  31 -DEG  HALF -ANGLE  SHARP 

CONE .  29 

10.  SUMMARY  OF  PROPOSED  CORRELATION 
PARAMETERS  UNDER  ENTROPY-LAYER¬ 
SWALLOWING  CONDITIONS .  37 

NOMENCLATURE .  46 


5 


AEDC-TR-75-2 


1.0  INTRODUCTION 


As  has  been  pointed  out  by  Masek  and  Forney  (Ref.  1),  two  signifi¬ 
cant  problems  in  the  application  of  existing  boundary-layer  technology 
to  the  National  Aeronautics  and  Space  Administration  (NASA)  Space 
Transportation  System  Space  Shuttle  Orbiter  are  (1)  the  choice  of  a 
boundary -layer  transition  criterion  and  (2)  the  method  used  to  determine 
boundary -layer  edge  conditions.  The  way  in  which  these  problems  are 
solved  has  a  large  effect  on  the  resulting  surface  temperature,  heating 
rate,  and  heat  load  distributions  along  the  lower  (windward)  surface  of 
the  vehicle  under  atmospheric  reentry  conditions  (Ref.  2). 

The  present  work  addresses  the  second  of  the  above-noted  problems, 

i.e,,  determination  of  boundary-layer  edge  conditions  for  the  specific 
case  of  hypersonic  flow  over  a  blunt-nosed  body  where  account  must  be 
taken  of  entropy-layer  swallowing  by  the  boundary  layer.  As  is  well 
known,  hypersonic  wind  tunnel  data  must  be  appropriately  extrapolated 
to  flight  conditions  by  analytical  methods  if  they  are  to  be  effectively 
utilized  in  vehicle  design.  In  keeping  with  this,  the  primary  objective 
of  the  present  study  is  the  development  of  correlation  parameters  which 
enable  perfect  gas  hypersonic  wind  tunnel  data  to  be  properly  scaled  to 
equilibrium  real  gas  flight  conditions  with  respect  to  laminar  boundary- 
layer  edge  quantities,  including  effects  of  entropy-layer  swallowing,  A 
secondary  objective  is  the  development  of  similarity  parameters  for 
correlation  of  surface  heat  transfer,  including  effects  of  entropy-layer 
swallowing.  The  physical  application  of  the  present  study  concerns  the 
extrapolation  to  flight  of  the  30- deg  angle -of -attack  laminar  boundary- 
layer  edge  condition  and  surface  heat -transfer  data  obtained  during  the 
recent  NASA /Rockwell  International  OH9  and  OH4B  tests  conducted  in 
the  AEDC  von  Karman  Gas  Dynamics  Facility  (VKF)  Hypersonic  Wind 
Tunnel  (B)  on  the  139  Space  Shuttle  Orbiter  configuration.  The  approach 
used  in  the  present  work  may  be  summarized  as  follows: 

1.  A  realistic  flow -field  model  amenable  to  existing  AEDC/ 

VKF  analysis  techniques  was  selected. 

2.  The  laminar  boundary-layer  flow  (including  entropy- 
layer  swallowing)  was  calculated  for  a  range  of  perfect 
gas  wind  tunnel  and  equilibrium  real  gas  flight  conditions, 

3.  The  results  were  correlated  in  terms  of  aerodynamic 
parameters  that  can  be  effectively  used  in  design  studies. 
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4,  Chemical  nonequilibrium  effects  of  the  inviscid  flow 
field  were  estimated  in  order  to  assess  the  applica- 
bility  of  the  equilibrium  real  gas  flight  calculations. 


2.0  FLOW-FIELD  MODEL 


Calculations  from  the  flow-field  model  developed  in  Ref,  3  for  sharp- 
edged  delta  wings  at  incidence  in  hypersonic  flow  indicate  that  an  angle- 
of -attack  range  exists  where  the  windward  centerline  flow  field,  both 
inviscid  and  viscous,  is  well  approximated  by  a  sharp  cone  at  zero  inci¬ 
dence  having  a  half-angle  equal  to  the  angle  of  attack.  In  the  case  of  an 
80-deg-sweep  delta  wing  at  free-stream  Mach  number  9.6,  this  angle- 
of-attack  range  extends  from  approximately  15  to  33  deg,  as  shown  by 
Fig.  74  in  Ref.  3,  The  windward  surface  of  the  Rockwell  International 
139  Space  Shuttle  Orbiter  (see  Fig,  1)  at  values  of  x/L  from  0,1  to  0,6 

X 

L 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

I - 1 - \ - -1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

X 

L 

Figure  1.  Rockwell  International  139  Space  Shuttle  Orbiter  Configuration. 
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has  a  relatively  flat  bottom  and  a  sweep  angle  of  80  deg.  The  windward 
centerline  profile  has  a  nearly  continuous  monotonically  decreasing 
slope  from  the  vehicle  nose  to  an  x/L  value  of  0,  25.  From  this  point 
the  slope  becomes  constant  at  one  degree  (with  respect  to  the  vehicle 
reference  system)  up  to  an  x/L  value  of  0.8. 

Based  on  the  above  discussion,  a  reasonable  approximation  of  the 
windward  centerline  flow  field  for  the  139  Space  Shuttle  Orbiter  at  30 
deg  angle  of  attack  should  be  an  appropriate  axisymmetric  body  at  zero 
incidence.  A  logical  starting  point  in  the  search  for  a  matching  axi¬ 
symmetric  body  is  a  31-deg  asymptotic  half-angle  hyperboloid.  The 
continuous  geometry  of  the  hyperboloid  makes  it  ideal  for  numerical 
flow-field  calculations.  Use  of  a  trial-and-error  method  with  a  model 
scale  drawing  showed  that  a  31-deg  hyperboloid  with  a  nose  radius  to 
body  length  value  of  0.02635  has  an  excellent  windward  centerline 
profile  comparison  with  the  139  Space  Shuttle  Orbiter,  as  shown  in 
Fig.  2.  For  the  full-scale  vehicle  the  values  of  body  length  and  nose 
radius  are  1290.3  and  34.0  in.,  respectively;  for  a  0.0175-scale  wind 
tunnel  model,  these  values  become  22.5  8  and  0.595  in.,  respectively. 


X 


Figure  2.  Comparison  between  139  Orbiter  windward  centerline  surface 
profile  at  30  deg  angle  of  attack  and  a  31-deg  asymptotic 
half-angle  hyperboloid. 
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Note  that  the  coordinate  system  defined  in  Fig.  2  is  aligned  with  the 
free-stream  velocity  vector  and  that  the  origin  is  at  the  body  point  where 
the  surface  is  normal  to  the  free-stream  velocity.  Also  shown  in  Fig.  2 
is  the  vehicle  coordinate  x/L.  The  correspondence  of  these  coordinates 
and  the  surface  distance  coordinate  s/r-^  is  given  in  Table  1. 


Table  1.  Coordinate  Correspondence 


x/L 

x'/L 

s/r^ 

0.  003 

0 

0 

0.  01 

0. 0061 

0,  74 

0.  02 

0. 0147 

1.27 

0,  05 

0. 0407 

2.  59 

0.  10 

0,  0840 

4.  62 

0.  15 

0. 1273 

6,  59 

0.  20 

0.  1706 

8.  53 

0.  3 

0.  2572 

12.  40 

0.  4 

0.  3438 

16.25 

0,  5 

0. 4304 

20.  11 

0.  6 

0.  5170 

23.  95 

0.  7 

0. 6036 

27.  79 

0.  8 

0. 6902 

31. 63 

Note:  a  =  30  deg 


Now  that  a  matching  axisymmetric  body  at  zero  incidence  has  been 
found,  flow-field  calculations  can  be  made  and  compared  to  hypersonic 
wind  tunnel  experimental  data  to  validate  the  proposed  flow-field  model. 
The  methods  used  are  described  in  following  sections. 

3.0  FLOW-FIELD  CALCULATION  METHODS 

3.1  BOUNDARY  LAYER  AND  INVISCID  FLOW 

Thre  theory  and  numerical  scheme  used  in  obtaining  the  present 
boundary-layer  results  are  based  on  the  work  of  Patankar  and  Spalding 
(Ref.  4).  In  this  approach  the  classical  boundary-layer  equations,  with 
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the  addition  of  the  transverse  curvature  terms,  are  expressed  in  a 
normalized  von  Mises  coordinate  system  and  solved  by  a  marching, 
implicit  finite-difference  procedure.  Although  the  primary  use  of 
the  work  of  Patankar  and  Spalding  has  been  in  the  analysis  of  turbu¬ 
lent  boundary  layers,  the  present  application  is  concerned  exclusively 
with  laminar  boundary  layers.  Further  details  concerning  the  digital 
computer  code  and  its  application  are  given  in  a  later  paragraph. 

Since  an  item  of  primary  interest  in  the  present  work  is  the  effect 
of  entropy  “layer  swallowing  on  laminar  boundary -layer  parameters,  a 
discussion  of  the  swallowing  process  and  its  analytical  treatment  is  in 
order.  The  specification  of  the  conditions  along  the  outer  edge  of  the 
boundary  layer  is  reasonably  straightforward  for  bodies  such  as  flat 
plates  and  sharp  cones;  however,  for  the  case  of  a  blunt  body  the  prob¬ 
lem  is  more  complex.  The  simplest  and,  therefore,  the  most  commonly 
employed  method  is  to  first  obtain  a  pressure  distribution  over  the  body 
of  interest  from  an  inviscid  method-of-characteristics  solution  or  New¬ 
tonian  theory.  Free-stream  stagnation  conditions  behind  a  normal  shock 
are  then  determined,  and  the  local  flow  conditions  along  the  boundary- 
layer  edge  are  found  by  isentropically  expanding  from  these  stagnation 
conditions  to  the  known  local  pressure,  thus  treating  the  boundary-layer 
edge  as  an  isentropic  surface.  This  method  is  satisfactory  for  the  for¬ 
ward  portion  of  a  blunt  body;  however,  as  the  flow  proceeds  along  the 
body  and  the  boundary  layer  grows  because  of  entrainment  of  mass,  the 
high  entropy  portion  of  the  flow  which  crossed  the  essentially  normal 
portion  of  the  bow  shock  is  swallowed  by  the  boundary  layer.  The  flow 
along  the  edge  of  the  boundary  layer  on  the  aft  portions  faf  the  body  will 
then  have  passed  through  an  oblique  part  of  the  bow  shock  and  will  be 
in  a  different  state  than  had  it  passed  through  a  normal  shock  (see 
Fig.  3). 


Assuming  the  pressure  along  the  outer  edge  of  the  boundary  layer 
on  a  blunt  body  to  have  the  inviscid  surface  value*,  the  determination 
of  the  local  edge  flow  conditions  may  be  improved  by  taking  into  con¬ 
sideration  the  inclination  of  the  bow  shock  where  the  local  flow  stream¬ 
line  crossed  the  shock.  The  point  at  which  the  flow  along  the  edge  of 
the  boundary  layer  crossed  the  shock  can  be  determined  by  matching  the 
mass  flow  in  the  boundary  layer  at  a  given  location  to  the  free-stream 


*This  assumption  should  be  valid  as  long  as  the  pressure  gradient 
normal  to  the  body  surface  is  small. 
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mass  flow  in  a  cylinder  with  radius  extending  out  to  the  location  to  be 
determined.  Referring  again  to  Fig.  3,  this  may  be  expressed  as 


2ttv  pu  dy 


(1) 


Note  that  y^  is  the  thickness  of  the  region  over  which  the  boundary- 
layer  equations  are  being  solved.  After  y^  is  founds  the  shock  inclina¬ 
tion  at  that  point  can  be  determined  and  the  flow  conditions  along  the 
boundary  layer  at  the  corresponding  body  location  can  be  computed  by 
crossing  the  oblique  shock  at  yg  with  the  free-stream  flow  and  allowing 
that  flow  to  expand  isentropically  to  the  known  local  boundary -layer 
edge  pressure.  Naturally,  a  swallowing  analysis  such  as  described 
above  requires  that  the  shape  of  the  bow  shock  be  known  in  addition  to 
the  body  surface  pressure.  In  the  present  calculations,  these  data  for 
both  real  and  perfect  gas  conditions  were  obtained  from  a  method-of- 
characteristics  solution  using  the  digital  computer  code  developed  by 
Inouye,  Rakich,  and  Lomax  (Ref.  5),  For  the  perfect  gas  cases,  the 
finite-difference  method  of  Aungier  (Ref.  6)  was  used  to  determine  the 
blunt  nose  subsonic  region  solutions  for  input  to  the  downstream  perfect 
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gas  method“of "Characteristics  analysis »  For  the  real  gas  (equilibrium 
air)  cases,  the  inverse  method  of  Lomax  and  Inouye  (Ref.  7)  was  used 
to  determine  the  blunt  nose  subsonic  region  solutions  for  input  to  the 
downstream  real  gas  method-of-characteristics  analysis. 

When  entropy  swallowing  is  being  considered,  the  definition  of  the 
outer  edge  of  the  boundary  layer  must  be  reconsidered.  This  is  neces¬ 
sary  because  the  velocity  and  temperature  gradients  are  not  zero  at  the 
outer  edge  of  the  boundary  layer  but  have  values  associated  with  the 
inviscid  flow  field.  This  is  the  result  of  the  combination  of  the  entropy- 
layer  swallowing  and  the  use  of  the  nonsimilar  boundary-layer  equations. 
To  treat  this  situation,  one  may  define  the  boundary-layer  thickness  in 
terms  of  the  total  enthalpy,  since  the  gradient  of  this  quantity  does  go 
to  zero  at  the  boundary -layer  edge.  This  approach,  proposed  by  Levine 
(Ref.  8),  was  used  herein.  In  particular,  the  region  treated  by  the 
boundary -layer  equations  was  adjusted  so  that  y^  was  a  minimum  of 
1, 1  times  the  value  of  y  where 


=  0.999  (2) 


Consideration  of  the  total  enthalpy  profile  behavior  also  allows  a 
boundary -layer  thickness  to  be  defined  in  terms  of  the  deviation  of  H/H„ 
from  unity.  A  total  enthalpy  boundary-layer  thickness,  has  been 
defined  as  the  distance  from  the  body  surface  to  the  point  where  H/H^^ 
differs  from  unity  by  0,005,  when  approached  from  the  outer  part  of  the 
flow.  For  hot-wall  cases  with  overshoot  in  H/Hg,  (i.e,,  values  greater 
than  one),  this  means  that  is  at  =  1,005  where  dH/dy  <  0,  For 

cases  with  no  overshoot,  it  merely  means  that  Sjj  is  at  H/H^  =  0.995, 

To  solve  parabolic  partial  differential  equations  such  as  the 
boundary -layer  equations  treated  herein,  it  is  necessary  to  specify 
both  the  boundary  conditions  on  the  dependent  variables  and  the 
initial  values  of  the  dependent  variables  in  the  cross -stream  direc¬ 
tion,  Experience  has  shown  that  any  choice  of  initial  profiles  which 
is  not  entirely  unreasonable  can  be  used  with  a  negligible  effect  on 
the  downstream  results  if  the  computations  are  begun  near  the  leading 
edge  of  a  sharp  body  or  the  stagnation  point  of  a  blunt  body.  This 
experience  has  been  used  to  obtain  laminar  boundary-layer  solutions 
in  the  vicinity  of  the  stagnation  point  of  a  blunt  body  by  beginning  the 
solution  near  the  stagnation  point  with  assumed  profiles  and  using  a 
number  of  very  small  steps  in  the  streamwise  direction.  With  the 
outer-edge  boundary  condition  determined  by  the  streamline -swallowing 
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technique  given  above,  the  numerical  computation  proceeds  downstream 
using  the  marching,  implicit  finite -difference  method  of  solution  dis¬ 
cussed  previously.  Note  that  this  approach  properly  considers  the  non¬ 
similar  growth  of  the  boundary  layer  along  the  body.  In  the  present 
analysis,  no-slip  boundary  conditions  are  applied  at  the  body  surface 
with  respect  to  velocity  and  temperature. 

The  computer  program  used  in  the  present  work  is  written  in  FOR¬ 
TRAN  IV  for  use  on  an  IBM  370/155  digital  computer  and  is  essentially 
a  highly  modified  version  of  the  code  originally  devised  by  Patankar  and 
Spalding  (Ref.  4).  This  program,  formulated  by  Mayne  and  Dyer  (Ref. 
9),  embodies  a  number  of  significant  modifications  and  extensions  to 
the  basic  technique  of  Patankar  and  Spalding.  These  include  elimination 
of  the  Couette  flow  analysis  at  the  body  surface  and  the  so-called  slip- 
value  scheme  in  favor  of  applying  the  basic  finite -difference  scheme 
across  the  entire  boundary  layer.  The  introduction  of  a  variable  cross - 
stream  step  size  which  is  smaller  near  the  wall  than  at  the  outer  edge  of 
the  boundary  layer  yields  accurate  solutions  while  still  permitting  the 
computations  to  proceed  efficiently.  Because  of  the  implicit  nature  of 
the  finite -difference  scheme,  no  numerical  stability  problems  are  to  be 
expected  from  this  approach;  experience  has  indicated  this  to  be  true. 

Calculations  for  the  real  gas  flight  conditions  in  the  present  study 
treat  the  flow  as  dissociating  air  in  thermodynamic  equilibrium  using 
the  thermodynamic  and  transport  property  correlation  formulas  of 
Cohen  (Ref,  10),  Under  hypersonic  wind  tunnel  conditions,  the  flow  is 
taken  to  be  thermally  and  calorically  perfect  air  with  a  constant  Prandtl 
number  of  0,7  in  conjunction  with  Sutherland's  viscosity  law. 

A  more  complete  discussion  of  entropy-layer  swallowing  by  laminar 
boundary  layers,  including  comparisons  of  results  from  boundary-layer 
analyses  which  include  and  exclude  entropy  swallowing  effects  relative 
to  a  fully  viscous  shock  layer  method,  is  given  by  Mayne  and  Adams  in 
Ref.  11, 


3.2  CHEMICAL  NONEQUILIBRIUM  STREAMTUBE 

The  theory  and  numerical  scheme  used  in  obtaining  the  present  in- 
viscid  chemical  nonequilibrium  results  are  based  on  the  work  of  Adams 
(Ref.  12)  for  oblique  shock  relaxation  in  air.  For  a  given  shock  angle, 
chemically  frozen  real  gas  Rankine-Hugoniot  oblique  shock  relations 
are  used  to  determine  the  gas  state  immediately  behind  the  shock.  A 
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chemical  nonequilibrium  streamtube  program  based  on  the  implicit 
integration  technique  of  Lomax  and  Bailey  (Ref.  13)  is  used  to  integrate 
the  governing  streamtube  equations  along  a  constant  pressure  stream¬ 
line  allowing  chemical  kinetics  to  start  immediately  behind  the  shock 
and  continue  downstream,  A  schematic  of  the  oblique  shock  relaxation 
process  is  given  in  Fig.  4. 


Figure  4.  Chemical  Nonequilibrium  Oblique  Shock  Relaxation. 


The  basic  gas  model  is  a  multicomponent  mixture  of  chemically 
reacting  perfect  gases  made  up  of  seven  species  (O,  O2,  N,  N2,  NO, 
NO"^,  and  e~).  The  gas  is  taken  to  be  in  vibrational  equilibrium  but 
in  chemical  nonequilibrium,  as  it  is  controlled  by  the  chemical  kinetic 
model  of  dissociating-ionizing  air  used  by  Cressv/ell,  et  al.  (Ref.  14). 
The  seven  chemical  reactions  are  given  in  Table  2;  the  corresponding 
forward  and  backward  reaction  rates  are  listed  in  Table  3,  while 
Table  4  gives  the  third-body  efficiencies  for  reactions  one,  two,  and 
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Table  2.  Chemical  Reactions 


Note:  represents  any  nonreacting  collision  partner 

that  may  be  present;  data  concerning  these  so- 
called  third  bodies  are  given  in  Table  4. 


Table  3.  Reaction  Rates 


Notes;  (1)  k  =  exp  (-E/RT) 

(2)  T  in  °K,  R  =  1.986  cal/gm-mole  °K 

O 

(3)  Units  are  gm-moles,  sec,  cm  ,  cal,  “K 
Reference:  Cresswell,  et  al.  (Ref.  14) 

Table  4.  Third-Body  Efficiencies* 


i 

*^2 

Ng 

N 

O 

NO 

1 

5 

2 

2 

25 

2 

2 

2 

5 

3 

5 

2 

3 

2 

2 

5 

5 

2 

^Multiplicative  factors  for  reaction  rates  involving 
nonreacting  collision  partners,  Mj^ 
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three  in  Table  2.  Thermodynamic  properties  (enthalpies  and  specific 
heats)  for  the  individual  species  are  taken  from  Blottner  (Ref.  15), 
whereas  the  mixture  viscosity  is  obtained  from  Wilke's  formula  using 
individual  species  properties  from  Ref,  15. 


4.0  FREE-STREAM  CONDITIONS 


4.1  WIND  TUNNEL 

Table  5  gives  a  set  of  three  free-stream  conditions  in  Tunnel  B 
for  which  data  were  taken  during  recent  tests  of  the  139  Orbiter  (Ref. 
16).  The  first  condition  corresponds  to  the  OH9  test  where  the 


Table  5.  Tunnel  B  Conditions 


Case 

Moo 

Reoo/ft 

RSoo, 

T  °R 

^O,  oo  •  ^ 

Test 

IT 

7.  92 

681,  614 

33,  794 

1332 

OH9 

2T 

8.  00 

3,  706,  270 

183,  757 

1351 

OH4B 

3T 

7.  92 

541, 605 

26,  853 

1263 

OH4B 

Note:  Reoo  «  is  based  on  a  nose  radius  of  0.04958  ft. 


boundary -layer  edge  conditions  were  measured  via  flow-field  probe 
surveys.  Model  wall  temperatures  were  also  measured  at  this  test 
condition;  these  values  are  tabulated  in  Table  6,  The  other  two  free- 
stream  conditions  in  Table  5  correspond  to  the  maximum  and  minimum 
Reynolds  number  conditions  of  the  OH4B  heat -transfer  test.  It  should 
be  noted  that  these  are  the  maximum  and  minimum  Reynolds  numbers 
available  in  Tunnel  B,  Model  wall  temperatures  for  these  conditions 
are  tabulated  in  Table  5. 

Note  the  case  designation  IT,  2T,  and  3T  in  Tables  5  and  6.  This 
designation  will  be  used  in  the  remainder  of  this  report  to  identify  the 
wind  tunnel  conditions, 

4.2  FLIGHT 

Table  7  gives  a  set  of  ten  altitude -velocity  flight  conditions  as 
specified  by  Rockwell  International  personnel  based  upon  the  Rockwell 
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Table  6.  Model  Wall  Temperature  Ratio  Distributions  under 
Tunnel  B  Conditions 


Note:  Above  values  are  T^/Tq  with  from  Table  5. 


Table  7.  Rockwell  International  SSV  Trajectory  No.  14040 


Case 

Time, 

sec 

Altitude, 

ft 

Velocity, 
ft /sec 

Moo 

RCoo/ft 

rjsj 

1 

400 

269, 915 

26, 006 

29.  41 

2,  610 

7,  395 

2 

496 

252, 149 

25,  576 

28.  01 

6,  064 

17, 181 

3 

608 

243,  742 

24,  727 

26.  38 

8,  193 

23,  213 

4 

704 

239, 322 

23,  864 

25.  14 

9,  390 

26,  605 

5 

800 

233,  529 

22, 897 

23.  73 

11, 214 

31, 773 

6 

896 

225,  698 

21, 767 

22.  09 

14,  177 

40, 168 

7 

1008 

212, 699 

20,  119 

19.  74 

20,  525 

58,  154 

8 

1104 

195,  803 

18,  148 

17.  22 

32, 755 

92,  806 

9 

1200 

186, 076 

15,  665 

14.  70 

39, 921 

113, 109 

10 

1296 

178, 440 

13, 310 

12.  39 

44, 194 

125, 216 

Notes:  (1)  RBoo  rjvji®  based  on  a  nose  radius  of  2.  83333  ft. 

(2)  For  this  trajectory  29.5  deg  <  a  <  31.6  deg;  all  calculations 
assume  a  =  30  deg. 
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International  SSV  Trajectory  No,  14040,  Graphical  presentation  of  this 
trajectory  is  given  in  Fig_  5,  The  corresponding  free-stream  prop¬ 
erties  were  taken  from  the  altitude -velocity  tables  of  Ref,  17,  which 
are  based  on  the  196  2  U,  S,  Standard  Atmosphere,  Note  the  numerical 
case  designation  numbers  1  through  10  in  Table  7  and  Fig,  5,  This 
designation  will  be  used  in  the  remainder  of  this  report  to  identify  the 
flight  conditions. 


Case 


Velocity,  ft/sec  x  10"^ 

Figure  5.  Rockwell  International  SSV  Trajectory  No.  14040. 

Radiative  equilibrium  wall  temperatures  for  a  laminar  boundary 
layer  along  the  139  Orbiter  windward  centerline  as  specified  by  Rockwell 
International  personnel  based  upon  the  Rockwell  International  SSV  Tra¬ 
jectory  No,  14040  are  given  in  Table  8,  These  wall  temperature  dis¬ 
tributions  were  input  to  the  boundary -layer  calculations  to  properly  take 
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into  account  variable  wall  temperature  effects  on  the  entropy-layer 
swallowing  phenomenon  under  flight  conditions. 

Table  8.  Centerline  Radiative  Equilibrium  Wall  Temperature 

Distributions  for  Rockwell  International  SSV  Trajectory 
No.  14040 


— .v^_^x/L 

0.  005 

0,  01 

0.  02 

0.  05 

0.  1 

0,  2 

0.  4 

0.  6 

0.  8 

1 

2250 

2119 

1944 

1713 

1553 

1472 

1415 

1323 

1199 

2 

2533 

2392 

2199 

1944 

1782 

1726 

1661 

1553 

1412 

3 

2594 

2451 

2254 

1995 

1830 

1772 

1706 

1594 

1451 

4 

2587 

2445 

2250 

1992 

1827 

1759 

1694 

1583 

1441 

5 

2585 

2444 

2250 

1993 

1828 

1751 

1688 

1577 

1435 

6 

2589 

2449 

2256 

1998 

1833 

1745 

1682 

1572 

1430 

7 

2596 

2457 

2264 

2006 

1841 

1735 

1673 

1563 

1422 

8 

2590 

2454 

2263 

2007 

1843 

1718 

1657 

1547 

1408 

9 

2383 

2259 

2083 

1847 

1694 

1541 

1486 

1392 

1265 

10 

2116 

2006 

1848 

1636 

1477 

1334 

1286 

1203 

1091 

Note;  Wall  temperature  values  are  in  °F. 


5.0  FLOW-FIELD  MODEL  VERIFICATION  UNDER 
HYPERSONIC  WIND  TUNNEL  CONDITIONS 


In  order  to  verify  the  applicability  of  the  axisymmetric  flow-field 
model  postulated  in  Section  2.0,  comparisons  have  been  made  between 
the  present  analysis  techniques  and  experimental  measurements  taken 
on  the  139  Orbiter  under  the  Tunnel  B  conditions  of  Tables  5  and  6  as 
discussed  in  Section  4, 1.  Figure  6  presents  the  measured  and  calcu¬ 
lated  inviscid  flow-field  properties.  The  surface  pressure  distribution 
shown  in  Fig.  6  reveals  relatively  good  agreement  between  measured 
and  calculated  values.  The  same  is  true  of  the  shock  shape  comparison 
shown  in  Fig.  6.  Note  that  the  differences  between  measured  and  cal¬ 
culated  shock  shape  coordinates  are  consistent  with  the  differences  in 
body  shape  coordinates  shown  previously  relative  to  Fig,  2. 
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1.0 


Po 


o  Experimental  Data,  AEDC/VKF  Tunnel  B, 
139  Orbiter,  a  =  30  deg,  0H9  Test 


—  Inviscid  Calculation,  Case  IT  (See  Table 
31-deg  Hyperboloid 

"o  5  C3  o  -0—0 


5). 


0  0.2  0.4  0.6  0.8  1.0 

x/L 


Calculated  Shxk 


* 

L 


Figure  6.  Comparison  between  hyperboloid  surface  pressure  and 

shock  shape  calculations  and  AEDC/VKF  Tunnel  B  data. 


A  comparison  between  measured  and  calculated  (based  on  the  pres- 
ent  entropy-layer-swallowing  laminar  boundary-layer  analysis)  local 
boundary -layer  edge  parameters  is  given  in  Fig.  7,  The  boundary- 
layer  edge  is  defined  as  in  the  discussion  in  Section  3,1;  i.e.,  all  edge 
conditions  (both  experimental  and  calculated)  are  based  on  a  total 
enthalpy  definition  of  the  boundary-layer  thickness,  Each  of  the 

parameters  shown  in  Fig,  7,  namely  boundary -layer  thickness,  edge 
Mach  number,  and  edge  Reynolds  number,  reveals  good  agreement  be¬ 
tween  experiment  and  calculations.  It  is  important  to  observe  the 
entropy-layer-swallowing  phenomenon  as  typified  by  the  edge  Mach  num¬ 
ber  and  edge  Reynolds  number  distributions  relative  to  the  corresponding 


21 


AEDC-TR-75-2 


limiting  values  for  inviscid  conical  flow  and  inviscid  normal  shock- 
isentropic  expansion.  Entropy -layer  swallowing  is  of  importance 
over  the  entire  body  for  this  particular  configuration  and  laminar 
flow  condition. 


o  Experimental  Data,  AEDC/VKF  Tunnel  B, 

139  Orbiter,  a  =  30  deg,  0H9  Test 

■  Present  Calculation,  Case  IT 

X  Present  Calculation,  Inviscid  I sentropic 
Expansion  from  the  Stagnation  Point 

+  Present  Calculation,  Inviscid  Conical  Flow 


L 

Figure  7.  Comparison  between  hyperboloid  boundary-layer  edge 
parameter  calculations  and  AEDC/VKF  Tunnel  B  data. 
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Calculated  and  measured  boundary -layer  profiles  are  compared  in 
Fig.  8  for  the  body  location  x/L  =  0.5.  Indicated  on  each  figure  is  the 
calculated  boundary-layer  thickness,  Sjj,  based  on  the  total  enthalpy 
definition.  As  can  be  seen  from  Fig.  8a,  pitot  pressure  measurements 
alone  cannot  be  used  to  experimentally  define  the  viscous-inviscid  inter¬ 
face  (i.e.,  the  boundary-layer  thickness)  under  hypersonic  flow  condi¬ 
tions  where  entropy-layer  swallowing  is  important.  Total  temperature 
measurements,  such  as  are  shown  in  Fig,  8b,  are  required  for  such 
flows.  It  should  be  noted  that  the  experimental  measurements  shown  in 
Figs.  8a  and  b  have  been  corrected  for  pressure  lag  error  and  radia¬ 
tion  losses,  respectively,  as  described  in  Ref.  16. 

Local  Mach  number,  tangential  velocity,  and  static  temperature 
profiles  are  presented  in  Figs.  8c,  d,  and  e,  respectively.  None  of 
these  profiles  reveals  any  indication  of  boundary-layer  thickness 
under  entropy-layer-swallowing  conditions.  The  experimental  mea¬ 
surements  shown  in  these  figures  represent  calculated  values  based 


a.  Pitot  Pressure  b.  Total  Temperature 


Figure  8.  Comparison  between  hyperboloid  boundary-layer  profile 
calculations  and  AEDC/VKF  Tunnel  B  data. 
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Body  Location  x/L  -  0. 5 


M 


c.  Mach  Number 


Figure  8.  Continued. 


upon  the  corrected  pitot  pressure  and  total  temperature  measurements 
of  Figs.  8a  and  b  in  conjunction  with  the  assumption  of  constant  static 
pressure  normal  to  the  surface. 

In  general,  agreement  between  calculated  and  measured  profile 
quantities  in  Figs,  8a  through  e  is  reasonably  good.  This  agreement 
further  illustrates  the  validity  and  applicability  of  the  axisymmetric 
flow -field  model  under  the  present  hypersonic  wind  tunnel  conditions. 

This  agreement  is  a  direct  function  of  the  good  shock  shape  simulation 
by  the  hyperboloid  geometry  shown  previously  in  Fig.  6, 

The  measured  (via  the  thin- skin  thermocouple  technique)  and  cal¬ 
culated  surface  heat-transfer  distributions  as  reflected  through  h/hj,gf 
are  given  in  Fig,  9,  The  comparison  is  quite  good  except  in  the  region 
of  0,15  <  x/L  <  0.30  where  the  measured  heat-transfer  rates  are  slightly 
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Body  Location  x/L  =  0.5 

o  Experimental  Data,  AEDC/VKF  Tunnel  B, 
139  Orbiter,  a  =  30  deg,  OH9  Test 

Present  Calculation,  Case  IT 


e.  Static  Temperature 
Figure  8.  Concluded. 


lower  than  the  calculated  values.  This  is  consistent  with  the  surface 
pressure  distribution  shown  previously  in  Fig.  6. 

Note  that  in  the  lower  portion  of  Fig.  9  the  boundary -layer  edge 
Mach  number  and  edge  Reynolds  number  at  the  onset  of  boundary-layer 
transition  are  indicated  based  upon  the  present  entropy -layer-swallowing 
laminar  boundary-layer  calculations.  In  the  absence  of  experimental 
flow-field  data  at  a  particular  test  condition,  an  entropy -layer-swallow¬ 
ing  boundary-layer  calculation  is  the  only  way  in  which  these  values  can 
be  obtained. 

The  results  shown  in  Figs.  6,  7,  8,  and  9  clearly  indicate  that  appli¬ 
cation  of  the  axisymmetric  flow-field  model  developed  in  Section  2,0  and 
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the  entropy-layer-swallowing  boundary-layer  calculation  methods  pres¬ 
ented  in  Section  3, 1  results  in  realistic  values  of  measurable  flow-field 
and  surface  parameters  for  the  139  Orbiter  at  30  deg  angle  of  attack 
under  hypersonic  wind  tunnel  conditions.  It  is  therefore  reasonable  to 
assume  that,  by  using  equilibrium  air  properties  in  the  present  calcula¬ 
tion  methods,  realistic  flow-field  calculations  can  be  obtained  under 
real  gas  flight  conditions,  if  the  flow  is  indeed  close  to  equilibrium. 


0  0.2  0.4  0.6  0.8  1.0 


L 


Figure  9.  Comparison  between  hyperboloid  surface  heat-transfer 
calculations  and  AEDC/VKF  Tunnel  B  data. 
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6.0  CORRELATION  OF  BOUNDARY-LAYER  PARAMETERS  INCLUDING 
EFFECTS  OF  ENTROPY-LAYER  SWALLOWING 


As  stated  in  Section  1,0,  the  most  important  objective  of  the  pres^ 
ent  investigation  is  the  development  of  correlation  parameters  which 
enable  hypersonic  wind  tunnel  data  to  be  properly  scaled  to  flight  with 
respect  to  laminar  boundary -layer  edge  quantities  under  flow  conditions 
where  entropy-layer  swallowing  must  be  taken  into  account.  Of  second¬ 
ary  importance  is  the  development  of  laminar  boundary-layer  heat- 
transfer  parameters  which  properly  reflect  entropy-layer  swallowing 
effects.  In  both  cases  the  chosen  parameters  must  be  simple  to  use 
(in  an  engineering  sense)  and  accurately  account  for  real  gas  effects 
as  encountered  in  flight. 


Following  Rotta  and  Zakkay  (Ref.  18),  whose  work  is  based  on 
earlier  work  by  Ferri  (Ref.  19),  Zakkay  and  Krause  (Ref.  20),  and 
Rubin  (Ref,  21),  there  is  sound  theoretical  justification  for  choosing 
the  following  form  of  a  similarity  parameter  (denoted  hereafter  as  S) 
for  the  ent  ropy -layer -swallowing  process 


S 


s/r. 


Re 


1/3 


(3) 


where  s  is  surface  distance  from  the  stagnation  point  (see  Fig,  3),  rj^ 
is  the  nose  radius  of  the  body,  and  Re^^  ^  free-stream  Reynolds 

number  based  on  the  nose  radius,  rj^.  The  1/3  power  on  Re^^^  in  the 

denominator  of  Eq.  (3)  is  due  to  the  theoretical  analysis  of  Ref,  18  as 
mentioned  above.  As  noted  in  Ref,  18,  the  dependence  of  the  laminar 
boundary-layer  outer  edge  properties  on  the  nose  radius  vanishes  when 
these  properties  are  expressed  in  terms  of  the  similarity  parameter  S 
under  flow  conditions  where  entropy-layer  swallowing  is  importarit,  A 
similar  finding  holds  with  respect  to  laminar  boundary-layer  heat 
transfer,  as  shown  in  Ref,  19,  It  is  to  be  emphasized  that  Eq,  (3) 
represents  a  similarity  parameter  for  the  entropy -swallowing  process 
on  axisymmetric  bodies  of  revolution  at  zero  angle  of  attack  under 
either  real  or  perfect  gas  conditions,  as  can  be  seen  by  reference  to 
the  theoretical  analysis  on  pp.  508-509  of  Ref.  18,  The  similarity 
parameter,  S,  is  not  applicable  in  the  nose  region  of  the  body;  neither 
is  it  applicable  for  very  low  Reynolds  number  flows.  The  correspond¬ 
ence  of  physical  surface  distance,  s,  for  a  given  value  of  S  is  given  in 
Fig.  10  for  both  the  wind  tunnel  and  flight  conditions  (Tables  5  and  7/ 
respectively)  of  current  interest.  Because  of  the  above  statement  that 
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Case 

2T  10  7  IT  3T  3 


®«  ''N 

Figure  10.  Correspondence  of  physical  surface  distance,  s,  for  a 
given  value  of  the  similarity  parameter,  S. 

the  similarity  parameter,  S,  is  not  applicable  in  the  nose  region  of  the 
body  or  for  very  low  Reynolds  number  flows,  it  is  assumed  that  the 
range  of  current  interest  as  indicated  in  Fig,  10,  namely  0. 10  <  S  <  1,0, 
is  a  reasonable  choice  for  the  present  flow  conditions.  All  succeeding 
figures  which  employ  the  similarity  parameter,  S,  will  be  restricted  to 
this  range. 
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Examination  of  the  calculated  results  for  flight  cases  1  through  10 
and  wind  tunnel  cases  IT  through  3T  has  revealed  the  following  correla- 
tion  parameters  for  the  outer  edge  boundary-layer  parameters: 


Outer  Edge  Quantity 


Re^/ft 


Correlation  Parameter 

^e/^e,  NS 
^conical 
(Reg/ft)/(Ree^ 


Here,  the  subscript  NS  denotes  the  so-called  normal  shock  outer  edge 
condition  which,  for  a  given  body  location,  is  obtained  by  isentropically 
expanding  the  body  surface  inviscid  flow  from  the  stagnation  point  to 
the  local  surface  pressure.  The  conical  Mach  number,  M^onicab 
based  upon  the  inviscid  sharp  cone  surface  Mach  number  as  tabulated 
in  Table  9  for  both  real  gas  flight  conditions  and  perfect  gas  wind  tunnel 
conditions;  note  that  the  sharp  cone  is  a  31-deg  half-angle  type,  which 
corresponds  to  the  asymptotic  half-angle  of  the  hyperboloid.  Use  of  the 


Table  9.  Conical  Mach  Number  Based  on  Inviscid  Solution 
for  a  31-deg  Half-Angle  Sharp  Cone. 


Case 

Type 

Gas 

^conical 

1 

Real 

5.  202 

2 

5.  085 

3 

4.  993 

4 

4.  987 

5 

5.  051 

6 

5.  070 

7 

4.  919 

8 

4.  631 

9 

4.  220 

10 

3.  782 

1T-3T 

Perfect 

2.  887 
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so-called  normal  shock  outer  edge  conditions  in  a  normalized  sense 
represents  an  attempt  to  reduce  the  real  gas  effects  on  the  resulting 
parameter;  the  free-stream  Mach  number  and  the  conical  Mach  number 
represent  a  further  attempt  to  account  for  free-stream  Mach  number 
effects.  The  power  of  unity  on  the  conical  Mach  number  and  0.70  on 
the  free-stream  Mach  number  are  the  results  of  plotting  Mg /Mg 
[or  (Ree/ft)/(Reg  Ns/ft)^  versus  Mgonical  ^oo)  a-  common  value 
of  the  similarity  parameter,  S,  for  both  the  flight  and  wind  tunnel  cases. 

Logarithmic  plots  which  show  a  data  band  encompassing  flight  cases 

I  through  10  and  wind  tunnel  cases  IT  through  3T  are  presented  in  Figs. 

II  and  12  for  the  outer  edge  Mach  number  and  outer  edge  Reynolds  num¬ 
ber  correlation  parameters,  respectively.  A  mean  fit  to  the  data  band 


Figure  11.  Logarithmic  plot  of  the  laminar  boundary-layer  outer 
edge  Mach  number  correlation  parameter. 
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Figure  12.  Logarithmic  plot  of  the  laminar  boundary-layer  outer 
edge  unit  Reynolds  number  correlation  parameter. 


is  shown  on  each  figure  along  with  the  resulting  numerical  constants  for 
the  mean  fit.  Also  indicated  by  closed  circles  on  the  figures  are  the 
experimental  data  from  the  Tunnel  B  flow-field  survey  along  the  wind¬ 
ward  centerline  of  the  139  Orbiter.  It  is  important  to  observe  from 
Figs.  11  and  12  that  the  maximum  deviation  from  the  mean  fit  to  the 
data  band  is  only  approximately  ±10  percent,  which  clearly  reveals  the 
applicability  of  the  presently  proposed  correlation  parameters  under 
both  real  gas  flight  conditions  and  perfect  gas  wind  tunnel  conditions. 
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Under  hypersonic  flow  conditions  where  entropy -layer  swallowing 
effects  are  important,  classical  laminar  boundary-layer  parameters 
such  as  displacement  thickness  and  momentum  thickness  become  ill- 
defined  due  to  gradients  in  the  boundary-layer  velocity  and  static  den¬ 
sity  profiles  at  the  outer  edge  of  the  boundary  layer.  A  boundary -layer 
parameter  which  has  physical  meaning  under  these  conditions  is  a 
boundary-layer  thickness  based  upon  the  total  enthalpy  definition  of  the 
boundary -layer  edge,  as  discussed  previously  in  Section  3.1.  Classical 
laminar  boundary -layer  theory  for  a  sharp  cone  at  zero  angle  of  attack 
in  a  supersonic  or  hypersonic  flow  indicates  that  the  boundary-layer 
thickness,  6,  is  given  by  an  expression  of  the  form 

(S/s')  xTRee.s'  =  f(M^,etc.)  (4) 


where  s'  is  surface  distance  from  the  apex  of  the  cone.  In  an  analo¬ 
gous  manner,  it  is  reasonable  to  postulate  that  the  laminar  boundary- 
layer  thickness,  6h>  based  on  a  total  enthalpy  definition  under  condi¬ 
tions  where  entropy-layer  swallowing  is  important  can  be  correlated 
by  an  expression  of  the  form 


(Sjj/s)  \/Ree,s  =  etc.) 


(5) 


which  can  be  written  as 


(Sh  / 


Re 


1/6 


=  f(IVl  etc.) 


s/r. 


Re 


1/3 


=  f(M^,  etc.)  VS 


(6) 


where  is  based  upon  the  local  value  of  Reg /ft  corresponding  to 

the  local  value  of  the  similarity  parameter,  S.  A  logarithmic  plot  of 
the  above-proposed  correlation  parameter  for  total  enthalpy -defined 
boundary -layer  thickness  is  given  in  Pig.  13  as  a  function  of  all  ten 
real  gas  flight  conditions  and  the  three  perfect  gas  wind  tunnel  condi¬ 
tions.  A  weak  free-stream  Mach  number  influence  on  the  correlation 
parameter  is  shown  in  Fig.  13,  namely 


(Sh/i'n)  VR®e,rj^ 
„  1/6  .,0.058 
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Data  Band  for  Flight  Cases  1  through  10 
and  Wind  Tunnel  Cases  IT  through  3T 


Mean  Fit  to  Data  Band 

(^h/^n)  \/«%,r. 


‘N 


2.615  S 


0.44 


j^l/6  „0.058 


s/rj^ 


Figure  13.  Logarithmic  plot  of  the  laminar  boundary-layer  thickness 
correlation  parameter. 


which  is  based  upon  examination  of  the  free-stream  Mach  number  effect 
at  a  common  value  of  the  similarity  parameter,  S.  Examination  of  Fig. 
13  reveals  that  the  above-proposed  correlation  parameter  for  total 
enthalpy -defined  boundary-layer  thickness  yields  correlated  results 
within  ±10  percent  maximum  deviation  from  the  mean  fit  to  the  data 
band  for  both  real  gas  flight  conditions  and  perfect  gas  wind  tunnel  con¬ 
ditions.  Note  also  that  Fig.  13  contains  the  experimental  data  (indicated 
by  closed  circles)  from  the  Tunnel  B  flow-field  survey  along  the  wind¬ 
ward  centerline  of  the  139  Orbiter. 


Let  us  now  turn  our  attention  to  development  of  a  surface  heat- 
transfer  correlation  parameter  under  flow  conditions  where  entropy- 
layer  swallowing  effects  are  important.  Classical  laminar  boundary- 
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layer  theory  for  a  sharp  cone  at  zero  angle  of  attack  in  a  hypersonic 
flow  under  cold  wall  conditions  indicates  that  the  surface  heat  trans¬ 
fer  is  given  by  an  expression  of  the  form 

Ste  \/R®e,s '  =  f(M^,  etc.)  (8) 

where  s'  is  surface  distance  from  the  apex  of  the  sharp  cone  and  Stg  is 
a  Stanton  number  based  on  local  edge  and  wall  conditions;  i.e,, 

Ste  =  - ^ -  (9) 

In  an  analogous  manner,  it  is  reasonable  to  postulate  that  the  laminar 
boundary -layer  surface  heat  transfer  under  conditions  where  entropy - 
layer  swallowing  is  important  can  be  correlated  by  an  expression  of 
the  form 


Ste  etc.) 


CIO) 


which  can  be  written  as 


Ste 


e.r„ 


Re 


1/6 


f(M^,  etc.) 


s/r 


N 


Re 


1/3 


1(M^,  etc.) 


(11) 


where  Reg  is  based  on  the  local  value  of  RCg/ft  corresponding  to 

the  local  value  of  the  similarity  parameter,  S.  In  the  evaluation  of  the 
Stanton  number 


Ste 


(12) 


the  value  of  the  local  edge  mass  flux,  PgUe^  i^^st  be  determined  in 
a  manner  consistent  with  the  local  value  of  Rog/ft;  in  addition,  the 
local  value  of  the  wall  enthalpy,  h.^,  corresponding  to  the  local  value 
of  the  similarity  parameter,  S,  must  be  used  for  nonuniform  wall 
enthalpy  conditions. 


Determination  of  the  local  edge  mass  flux  PgUg  in  a  manner  con¬ 
sistent  with  the  local  value  of  Reg /ft  is  a  subtle  point,  especially  for 
real  gas  flow  conditions.  The  following  iterative  calculation  procedure 
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is  suggested  to  determine  a  consistent  value  of  PgU^  under  real  gas 
conditions: 

1. 


2. 


3. 


4. 


Under  perfect  gas  wind  tunnel  conditions,  exactly  the  same  calculation 
procedure  is  used  except  that  in  step  (3)  the  perfect  gas  equation  of  state 
is  used  to  determine  and  Sutherland's  viscosity  law  is  used  to  deter¬ 
mine  u  . 
e 

A  logarithmic  plot  of  the  above  proposed  correlation  parameter  for 
surface  heat  transfer  is  given  in  Fig.  14  as  a  function  of  all  ten  real  gas 
flight  conditions  and  the  three  perfect  gas  wind  tunnel  conditions.  It  is 
to  be  noted  that  no  free-stream  Mach  number  effect  could  be  discerned 
from  examination  of  the  data  in  the  manner  discussed  previously.  As 
can  be  seen  from  Fig.  14,  the  above  proposed  correlation  parameter  for 
surface  heat  transfer  yields  correlated  results  within  ±5  percent  maxi¬ 
mum  deviation  from  the  mean  fit  to  the  data  band  for  both  real  gas  flight 
conditions  and  perfect  gas  wind  tunnel  conditions.  The  key  points  in  this 
successful  application  are  the  determination  of  the  local  mass  flux  PgUg 
per  the  above  discussion  and  the  use  of  the  local  wall  enthalpy  in  the 
Stanton  number  expression. 


Find  Reg/ft  =  p^  Ug/Pg  at  a  given  value  of 
S  =  (s/r]^)/Rei^' from  the  correlation  curve. 

Guess  a  value  for  Ug  (a  reasonable  first  guess  is 
Ug  =  Uoo  cos  a)  and  calculate  h^  from  constancy  of 
total  enthalpy,  namely  hg  =  -  1/2  . 

Note  that  for  a  real  gas  (equilibrium  air),  the  Cohen 
(Ref.  10)  correlation  functions  give  pg  =  fi(p^,  hg)  and 
PgPg  =  ^2^Pw^  where  the  local  value  of  the  surface 
pressure  p^  is  used  consistent  with  the  value  of  S. 


With  U^,  p^,  „„„  ,g..g 


and  PgPp  determined  by  the  above 
procedure, '"compute 

Pe  Ug  p|  Ug 

Re^/ft  = 


P 

'^e  e 


and  compare  the  resulting  value  with  the  value  from 
step  (1)  above.  Iterate  on  value  of  Ug  using  steps  (2), 
(3),  and  (4)  until  the  procedure  converges  to  consistent 
values  of  ROg/ft  from  steps  (1)  and  (4). 
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Data  Band  for  Flight  Cases  1  through  10 
and  Wind  Tunnel  Cases  IT  through  3t 


Mean  Fit  to  Data  Band 


Figure  14.  Logarithmic  plot  of  the  laminar  boundary-layer  surface 
heat-transfer  correlation  parameter. 


A  summary  of  the  presently  proposed  laminar  boundary-layer 
correlation  parameters  under  entropy -swallowing  conditions  is  given 
in  Table  10.  As  shown  above,  these  parameters  are  applicable  under 
real  gas  flight  conditions  as  well  as  perfect  gas  wind  tunnel  conditions. 
Through  use  of  these  parameters,  laminar  boundary-layer  probe 
survey  data  and  surface  heat -transfer  data  taken  in  a  hypersonic  wind 
tunnel  can  be  properly  scaled  to  flight  conditions  with  respect  to 
boundary-layer  outer  edge  parameters,  boundary-layer  thickness,  and 
surface  heat  transfer  under  flow  conditions  where  entropy-layer  swal¬ 
lowing  is  important. 


7.0  ESTIMATE  OF  CHEMICAL  NONEQUILIBRIUM 
EFFECTS  ON  OUTER  EDGE  PARAMETERS 


All  of  the  foregoing  analysis  relative  to  flight  conditions  is  based 
upon  the  assumption  of  thermochemical  equilibrium  in  both  the  inviscid 


36 


Table  10.  Summary  of  Proposed  Correlation  Parameters 
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Abscissa 


Ordinate 


S  = 


s/r 


N 


Re 


1/3 


^e/^e,  NS 
^conical 

Reg/ft/Reg,  Ns/^t 


N 


Re 


1/6 

»,rN 


M, 


0.  058 


Re, 


1/6 


■N 


and  viscous  flow  fields.  In  order  to  examine  the  effects  of  chemical 
nonequilibrium  on  boundary-layer  outer  edge  parameters  under  flight 
conditions,  a  series  of  calculations  has  been  performed  for  the  case 
of  oblique  shock  relaxation  in  air  using  the  chemical  nonequilibrium 
streamtube  analysis  described  in  Section  3.  2  of  the  present  report. 
Three  sets  of  altitude -velocity  conditions  have  been  considered,  namely 


Altitude,  ft 


Velocity,  ft /sec 


250, 000 
220, 000 
190, 000 


25,  600 
21, 300 
16,  700 


which  are  reasonably  close  to  the  Rockwell  International  SSV  Trajectory 
No.  14040  given  in  Table  7  and  shown  in  Fig.  5.  At  each  of  these 
altitude -velocity  conditions,  three  shock  angles  have  been  considered, 
namely  ^s^ock  "  where  Oghock  defined  as  in  Fig. 

4..  This  choice  of  shock  angles  is  designed  to  simulate  inviscid  stream¬ 
line  swallowing  under  chemical  nonequilibrium  conditions.  For  example. 
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0ghock  ■  represents  a  fully  swallowed  streamline  on  a  body  at 

approximately  31  or  32  deg  angle  of  attack;  ^shock  ~  represents 

a  typical  streamline  in  the  region  where  swallowing  phenomena  are 
important. 

Summary  plots  of  the  local  Mach  number  and  local  unit  Reynolds 
number  are  given  in  Figs,  15  and  16  for  the  three  altitude-velocity 
conditions  of  present  interest.  Also  shown  in  these  figures  are  equi¬ 
librium  air  results  for  the  same  shock  angles  and  flow  conditions.  It 
is  to  be  noted  that  the  local  Mach  number  is  defined  in  terms  of  the  so- 
called  frozen  speed  of  sound  under  chemical  nonequilibrium  conditions 
and  the  equilibrium  speed  of  sound  under  chemical  equilibrium  condi¬ 
tions;  see  Section  5  in  Chapter  VIII  of  Vincent!  and  Kruger  (Ref,  22) 
for  an  excellent  discussion  of  frozen  and  equilibrium  speeds  of  sound. 
As  can  be  seen  from  Fig.  15,  chemical  nonequilibrium  effects  on  the 
local  Mach  number  are  almost  insignificant,  especially  for  the  higher 
shock  angles.  A  more  significant  chemical  nonequilibrium  effect  is 
shown  in  Fig,  16  with  respect  to  the  local  unit  Reynolds  number.  For 
the  35 -deg  shock  angle  condition,  the  equilibrium  value  of  the  local  unit 


Oblique  Shock  Relaxation  '  ' 

250, 000  ft  at  25, 600  ftisec 
Nonequilibrium  Streamtube  Calculation 
Seven  Species  Air  (0,  O2,  N,  N2,  NO,  NO'*’,  e  ) 


®shock' 


Equilibrium 

Air 


Distance  Behind  Shock,  ft 


a.  250,000  ft  at  25,600  ft/sec 

Figure  15.  Chemical  nonequilibrium  oblique  shock  relaxation  effects 
on  the  local  Mach  number. 
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Oblique  Shock  Relaxation 
220,000  ft  at  2 1,300  ft/sec 
Nonequilibrium  Streamtube  Calculation 
Seven  Species  Air  (0,  O2,  N,  N2,  NO,  NO"'',  e'l 


®shock' 


jJ _ I _ I _ I  I  I  I  I  I  I _ I _ I _ I  I  I  I  M  I _ I _ I . . 1  I 

0. 1  1. 0  10. 0  100. 0 
Distance  Behind  Shock,  ft 

b.  220,000  ft  at  21,300  ft/sec 


Oblique  Shock  Relaxation 
190, 000  ft  at  16, 700  ft/sec 
Nonequilibrium  Streamtube  Calculation 
Seven  Species  Air  (0,  O2,  N,  N2,  NO,  NO'*',  e') 


®shock' 


Equilibrium 

Air 


1 1  I  I  I  I  I  1 1 1 1  I  I  I  I  I  1 1 1 1  I  I  I  I  I  1 1 1 1 
0.1  1.0  10.0  100.0 
Distance  Behind  Shock,  ft 

c.  190,000  ft  at  16,700  ft/sec 
Figure  15.  Concluded. 
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Equilibrium 

Air 


a.  250,000  ft  at  25,600  ft/sec 


Oblique  Shock  Relaxation 
220, 000  ft  at  21, 300  ft/sec 


Distance  Behind  Shock,  ft 


Equilibrium 

Air 


b.  220,000  ft  at  21,300  ft/sec 

Figure  16.  Chemical  nonequilibrium  oblique  shock  relaxation  effects 
on  the  local  unit  Reynolds  number. 
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Oblique  Shock  Relaxation 
190, 000  ft  at  16, 700  ft/sec 


Equilibrium 

Air 


c.  190,000  ft  at  16,700  ft/sec 
Figure  16.  Concluded. 

Reynolds  number  is  approximately  40  percent  greater  than  the  chemical 
nonequilibrium  value  at  a  distance  of  100  ft  behind  the  shock;  this  state¬ 
ment  applies  to  all  three  altitude -velocity  conditions.  At  the  55-deg 
shock  angle  condition,  equilibrium  is  essentially  attained  at  100  ft  be¬ 
hind  the  shock  for  the  two  lowest  altitude  conditions. 

The  above  findings  indicate  that  chemical  nonequilibrium  effects 
will  tend  to  reduce  both  the  local  Mach  number  and  local  unit  Reynolds 
number  below  the  corresponding  equilibrium  air  values  for  the  same 
altitude -velocity  and  shock  angle  conditions.  The  largest  chemical 
nonequilibrium  influence  is  on  the  local  unit  Reynolds  number  at  the 
smaller  values  of  the  shock  angle.  Much  more  detailed  calculations 
using  chemical  nonequilibrium  bluntbody  flow-field  analyses  are  re¬ 
quired  to  fully  assess  nonequilibrium  effects  under  Space  Shuttle  flight 
conditions.  Until  such  knowledge  becomes  available,  equilibrium  air 
results  appear,  in  light  of  the  above  findings,  to  yield  reasonable 
estimates  under  flow  conditions  where  entropy-layer  swallowing  is 
important. 
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8.0  IMPLICATIONS  RELATIVE  TO  BOUNDARY-LAYER  TRANSITION 


The  present  work  emphasizes  the  development  of  correlation  pa¬ 
rameters  which  enable  hypersonic  wind  tunnel  data  to  be  properly 
scaled  to  flight  with  respect  to  laminar  boundary-layer  edge  quantities 
under  flow  conditions  where  entropy-layer  swallowing  must  be  taken 
into  account.  The  importance  of  these  correlation  parameters  is 
directly  linked  to  the  common  aerodynamic  practice  of  correlating 
laminar  boundary-layer  transition  onset  in  terms  of  local  boundary- 
layer  edge  properties.  One  example  of  such  a  correlation,  taken  from 
Fig,  45  in  Ref.  23,  gives 


^  0  09^ 

Ree,st  =  0.984  X 


where 


Reg, St  =  {Reg /ft)  St 


(13) 


(14) 


with  s^  being  the  surface  distance  from  the  nose  to  onset  of  boundary- 
layer  transition.  Needless  to  say,  accurate  estimates  of  Reg/ft  and 
Mg  are  required  in  order  to  apply  the  above  correlation  (which  is  based 
on  perfect  gas  hypersonic  wind  tunnel  data  for  various  Phase  B  Shuttle 
configurations)  under  real  gas  flight  conditions  where  entropy-layer 
swallowing  is  important.  The  presently  proposed  correlation  param¬ 
eters  satisfy  this  requirement  in  all  respects. 


The  correlations  developed  in  this  investigation  allow  boundary- 
layer  edge  conditions  and  transition  estimates  to  be  determined  for 
perfect  gas  wind  tunnel  conditions  and  equilibrium  real  gas  flight  con¬ 
ditions.  The  above  correlation  for  onset  of  boundary -layer  transition 
can  be  used  in  conjunction  with  the  present  oblique  shock  relaxation 
calculations  to  estimate  the  effects  of  chemical  nonequilibrium  on 
boundary -layer  transition  onset  distance,  s^.  Denoting  chemical  non¬ 
equilibrium  by  NEQ  and  chemical  equilibrium  by  EQ,  the  above  cor¬ 
relation  yields  the  ratio 


(st)EQ 


Me 


0.923' 


L  Reg/ft 


'NEQ 


0,923 


Reg/ft 


EQ 


(15) 
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under  the  assumption  that  the  correlation  (which,  recall,  is  based  on 
perfect  gas  hypersonic  wind  tunnel  results)  is  applicable  under  both 
chemical  equilibrium  and  nonequilibrium  conditions.  Evaluating  the 
chemical  nonequilibrium  parameters  at  a  distance  of  50  ft  behind  the 
shock  for  the  altitude -velocity  condition  of  220,  000  ft  and  21,  300  ft /sec 
yields  the  following: 


^shock* 


35 

45 

55 


1. 30 
1.10 
1.04 


These  results  indicate  that,  for  a  given  altitude -velocity  and  shock  angle 
condition,  the  effects  of  chemical  nonequilibrium  will  be  to  increase 
the  transition  onset  distance.  It  is  to  be  emphasized  that  this  finding  is 
based  on  one  particular  transition  correlation;  other  choices  for  the 
transition  correlation  may  yield  different  trends  with  respect  to  chemical 
nonequilibrium  effects  on  transition  onset  distance.  More  work  is  defi¬ 
nitely  needed  in  this  area.  The  recent  study  by  Lordi,  Vidal,  and  John¬ 
son  (Ref.  24)  represents  one  such  effort  similar  in  many  respects  to  the 
present  investigation. 


9.0  CONCLUDING  SUMMARY 


The  present  investigation  has  resulted  in  development  of  correla¬ 
tion  parameters  for  laminar  boundary-layer  edge  quantities  and  surface 
heat  transfer  along  the  windward  symmetry  plane  of  the  Rockwell  Inter¬ 
national  139  Space  Shuttle  Orbiter  configuration  at  30  deg  angle  of 
attack.  These  correlation  parameters  enable  hypersonic  wind  tunnel 
data  to  be  properly  scaled  to  flight  under  conditions  where  entropy- 
layer  swallowing  by  the  boundary  layer  must  be  considered.  The  pro¬ 
posed  correlation  parameters  are  simple  to  use  (in  an  engineering 
sense)  and  accurately  account  for  real  gas  (equilibrium  air)  effects  as 
encountered  in  flight.  Chemical  nonequilibrium  effects  have  been 
briefly  examined  via  oblique  shock  relaxation  in  order  to  ascertain  the 
validity  of  the  equilibrium  air  flow-field  model  used  in  the  present 
work.  Some  implications  of  the  proposed  correlation  parameters  on 
boundary-layer  transition  have  been  discussed. 
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NOMENCLATURE 


H 

Ho 

h 

he 

h 


w 
h 

^ref 

L 


M 


M, 


conical 
Me 

Me,  NS 


Total  enthalpy 
Free-stream  total  enthalpy 
Static  enthalpy 

Static  enthalpy  at  edge  of  boundary  layer 
Wall  static  enthalpy 

Heat  “transfer  coefficient,  -q^/ (Tq^  ^  ~  T^) 

Reference  heat-transfer  coefficient  based  on  Fay-Riddell 
correlation  and  a  scaled  1-ft  spherical  nose  radius 

Reference  length  for  139  Orbiter,  1290,  3  in.  full  scale, 
22.  58  in.  for  0.  0175-scale  model 

Local  Mach  number 

Inviscid  conical  surface  Mach  number 

Local  Mach  number  at  edge  of  boundary  layer 

Local  Mach  number  at  outer  edge  of  boundary  layer  based 
on  isentropic  expansion  from  the  stagnation  point 

Free-stream  Mach  number 
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r 


Free-stream  pitot  pressure 
Local  pitot  pressure 
Wall  static  pressure 
Wall  heat  “transfer  rate 
Radius  as  defined  in  Fig,  3 


r^^  Nose  radius  as  defined  in  Fig,  3 

Re/ft  Local  unit  Reynolds  number 


Reg /ft  Local  unit  Reynolds  number  at  edge  of  boundary  layer, 

^e^e  ^  ^e 


Ree  yfs  Local  unit  Reynolds  number  at  edge  of  boundary  layer  based 
on  isentropic  expansion  from  the  stagnation  point 


Re, 


Re 


e,  s 


Re 


e,st 


RCoo/ft 


Local  Reynolds  number  at  edge  of  boundary  layer  based  on 
nose  radius,  PeUg^p^/Pg 

Local  Reynolds  number  at  edge  of  boundary  layer  based  on 
surface  distance  from  sharp  cone  apex,  p^U^s'/Pg 

Local  Reynolds  number  at  edge  of  boundary  layer  based  on 
surface  distance  from  the  stagnation  point,  p^UgS/p^ 

Local  Reynolds  number  at  edge  of  boundary  layer  based  on 
surface  distance  from  the  stagnation  point  to  onset  of 
boundary-layer  transition,  PgUgS^/p^ 

Free-stream  unit  Reynolds  number,  p^'Ua^/^a> 


Free-stream  Reynolds  number  based  on  nose  radius, 
PcoU^r^/Poo 


S 

s' 


Entropy-layer-swallowing  similarity  parameter. 
Surface  distance  from  the  sharp  cone  apex 
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